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ABSTRACT 

We present observations of the Class protostar L1157-mm using the Combined 
Array for Research in Millimeter-wave Astronomy (CARMA) in 3 mm dust continuum 
and N2H+ line emission. In the N2H"'" line, we detect a large-scale envelope extended 
over a linear size of ~20,000AU flattened in the direction perpendicular to the outflow. 
This N2H"'" feature coincides with the outer envelope seen in the 8 /xm extinction by 
Looney et al. Meanwhile, the dust continuum traces the compact, nearly spherical 
structure of the inner envelope, where N2H"'" becomes depleted. This highly flattened 
N2H"^ envelope also shows dynamical signatures consistent with gravitational infall in 
the inner region, but a slow, solid-body rotation at large scales. This flattened structure 
is not a rotationally supported circumstellar disk; instead, it resembles a prestellar 
core both morphologically and kinematically, representing the early phase of a Class 
system. In this paper, we construct a simple model to interpret both the dust continuum 
and N2H+ emission and suggest a possible dynamical scenario for the overall properties 
of the envelope. 

Subject headings: stars: formation — techniques: interferometric — ISM: molecules — 
stars: individual (L1157) — 



1. Introduction 

Significant progress has been made both theoretically and observationally toward understand- 
ing lo w-mass star formation in the last few decades, and an evolutio nary sequence has been posited 



(e.g.. lAndre et al.ll200d : iMcKee Qstrikeill2007l : Evans et al.ll2009l ). The youngest protostars, the 



so-called Class sources, form from the preceding stage of a prestellar core and are deeply embed- 
ded in their natal collapsing envelope. As a link from prestellar cores to young protostars, these 
envelopes contain valuable physical clues to the initial conditions of the collapse process. While 
single-dish telescope observations alone are usually limited by large beam size, interferometry has 
allowed high resolution studies of the Class envelope. In particular, interferometric observa- 
tions combined with dust modeling using theoretically predicted structure can probe the physical 
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properties of the inner envelope (e.g.. iLoonev et al.ll2003l : IChiang et al.ll2008l ). Unfortunately, dust 
emission from the outer envelope is usually difficult to detect due to interferometer insensitivity to 
the large-scale structure or the low-surface brightness of the extended emission. Molecular lines, 
on the other hand, trace specific components of the envelope and can detect regions of lower den- 
sity. Further, as the line emission is dispersed over many velocity channels, the structures are 
less extended than the equivalent dust emission so less affected by interferometer resolving-out 
issues. However, uncertainty in abundances and chemistry in the envelope can make molecular line 
interpretation difficult. 



Since the first detection in I Turned (|1974| ) and confirmation in iThaddeus fc Turneii (|l975l ) 
N2H"'" has been a well-knowii inter stellar molecule and an excellent tracer of dens e cores (e.g. 



Caselli et al. 



2002 



Chen et al.ll2007l ). With a critical density of around 10^ cm ^ (e.g.. lDaniel et al 



20061 ). N2H"'" is sensitive to the outer envelope around Class young stellar objects (YSOs). Also, 
N2H"'" has a low depletion rate and remains in t he gas phase when other molecules such as CO 
have depleted onto grains (iBergin &: LangedllQQTI ). The reason for this was thought to be the low 
binding energy of its parent molecule N2, but r ecent experiments have shown that the binding 
energy of N2 is not much lower than that of CO (lOberg et al.ll2005l ). In fact, most of the nitrogen 
is in the atomic form, which has a binding energy lower than the molecular form, and therefore 
making N2H"'" harder to deplete (jMaret et al.l 120061 ). In this paper, we use N2II+ as an envelope 
tracer to reveal the gas in the outer protostellar envelope, as well as pursue a comparison between 
gas and dust. 

While most Class protostellar envelopes have rather complex structures, L1157-mm (also 



know n as L1157-IRS or IRAS 20386+6751) in the Cepheus fiare fe.g.. iKun et al.ll2008l : iKirk et al 
2003) has a highly fiattened and relati vely symmetric envelope 2' across seen by 8 absorption 



(^Loonev et al. 



2007 



Tobin et al 



2003 ) . The protostar is in the Class stage with an estimated 



age of 8-150 kyr ([Froebrichl |2005| . based on bolometric temperature, bolometric luminosity, and 



submillimeter luminosity). Different components of the dus t emission, including a c e ntral compact 



core and an extended envelope, have been observed (e.g., iGueth et al.l 119971 . l2003l : iBeltran et al 



2004 ). Also, large-scale outfiow s have been detected by various molecules, suggesting an inclination 



angle of 80° (IGueth et a! 



clumps (jBachiller et al 



2001 



1996 ) and a kinematic age of 15 kyr derived from the oldest pair of outfiow 



Arce et al.ll2008l . and reference therein). Fur thermore, detection of 



methanol has suggest ed active accretion in the embedded circumstellar disk (jGoldsmith et al.lll999l : 
Velusamv et al.ll2002l ). With features of a typical Class young stellar object, L1157 seems to be 
a good example of isolated low-mass star formation. 

Unfortunately, the distance to L1157 is uncertain. The distance of 440 pc, based on spec- 
troscopy and ph otometry of the illuminating star HD200775 of NGC 7023, is commonly used (see 



Kun et al.l l2008l ) , while the distance of the neighbori ng L1147/L1158 comp lex was determined to 
be 325 pc by the extinction-distance relationship in IStraizys et al.l ( 19921 ). Later, the Cepheu s 



cloud was found to have three characteristic distances 200 pc, 300 pc, and 450 pc in Ik^ jiggi) 
Multiple layers of clouds make it difficult to know the exact distance to a specific region. For easier 
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comparison with the study of iLoonev et alj (|200?1 ). we adopt 250 pc in this paper 



We present interferometric observati ons of L1157-min using the Combined Array for Research 
in Milhmeter-wave Astronomy fCARMA: lBock et al.ll200^ )R The dust continuum at 3 mm and the 



N2H^ gas emission are observed and studied. The observational setup and the results are presented 
in §2 and §3. We estimate the column density of the N2H"'" emission (§4.1) and compare it with 
the previous detection of the flattened envelope at 8 nm absorption (§4.2). The N2H"'" abundance 
is examined in §4.3. We construct a simple model to interpret both the dust continuum and the 
N2H"'" gas emission in §4.4. In addition, the gas kinematics of the large-scale envelope is discussed 
in §4.5. An overall link between gas and dust at multiple scales in this system is concluded in §4.6, 
following by a summary in §5. 



2. Observations 

The observations were carried out with CARMA in October 2008 in E configuration and March 
2009 in D configuration. The system temperature ranged from 150 to 300K (single sideband, SSB) 
during source observations. D and E configurations are the two most compact configurations of 
CARMA and give angular resolutions of around 5.5" and 10" at 3 mm, respectively. Because of 
the short antenna separations (8-66 m) in E configuration, antennas can shadow each other when 
observing targets at low elevations. To minimize shadowing, the elevations were always above 
30° throughout the observations. Also we checked and verified that no shadowed data of the source 
were included in the analysis. 

CARMA is composed of nine 6.1m Berkeley Illinois Maryland Association (BIMA) antennas 
and six 10.4m Owens Valley Radio Observatory (OVRO) antennas [^. At our observed frequency, 
the FWHM of the primary beams are 122" and 70" for the 6m and 10m dishes, corresponding to 
30,500AU and 17,500AU, respectively. The sensitivity decreases significantly outside the FWHM. 
To observe the extended feature of L1157, we perfor med a five-pointing mosaic observation across 



the disk- like structure found in iLoonev et al.l (|2007l ). The phase center of the central pointing is 



set at the position of the protostar, a = 20''39"'06^26 and 6 = 68°02'15"8, determined by the dust 
continuum peak from other high resolution CARMA observations (Chiang et al. in preparation). 

The correlator was set so that the N2H+ JFiF=101-012 component at 93.1763 GHz was ob- 
served with a 2 MHz bandwidth, which provided a velocity resolution of 0.098 km s^^. Dust con- 
tinuum was observed simultaneously with two 500 MHz bands. During the observations, 1927+739 
was observed every 18 minutes and used as the phase calibrator. Also, the bandpass calibrators 
were 3C454.3, 3C84, 3C345, and 1642+689 for different tracks, while the flux calibrators were 
Uranus and MWC349. The absolute flux uncertainty is around 10%; hereafter we only consider 
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the statistical un certainty. All ob servational data were reduced and imaged using the MIRIAD 
software package ( Sault et al. 19951 ) [§. For the maps shown in this paper, natural weighting is used 
for Fourier transforming the visibilities into the image space. 



3. Results 

3.1. Dust Continuum 

Figure [T][a) shows the 3 mm dust continuum map. With a beam size of 7.3" x 6.5" (1825 
AUxl625 AU) in the combined D- and E-array CARMA data, the dust emission is compact and 
nearly spherical. We used the MIRIAD task IMFIT to fit the dust emission with a Gaussian, and the 
results are given in Table [TJ We detect the spherical protostellar envelope, but the embedded early 
circumstellar disk is not resolved. The dust continuum seen by CA RMA is more compact than that 



seen by single-dish observations (e.g., 850 /um and 1.3 mm maps in IShirlev et al.ll2000l : iGueth et al 



20031 ). implying that the large-scale extended emission is resolved out by the interferometer. A more 



careful modeling of the dust envelope is done in §4.4. 



3.2. N2H+ Gas 

Figures [TJl^b) and[2]show the integrated intensity (zeroth moment map) and the velocity channel 
maps of the N2H"'" isolated hyperfine component Ji^iF= 101-01 2. The star marks the position of the 
central protostar. We detect large-scale N2H"'" emission across an elongated region perpendicular 
to the outflow direction (see the gray-scale image in Figure ^a) for the outflow orientation) . The 
disk-like structure extends over 80", which corresponds to a linear size of ~20,000 AU assuming 
a distance of 250 pc. Due to the large size, it is unlikely to be a rotationally supported disk (see 
§4.5.2). Hence, we call it a flattened envelope in this paper. 

The N2H"'" gas shows interesting dynamics of the flattened envelope in FigureEl Systematically, 
the east wing has lower velocity than the west wing, suggesting rotation or bulk motion. The 
dynamics of the system will be discussed more in later sections. The spatial distribution is clumpy. 
In particular, the extended emission toward the east-south roughly follows the interface of the 
envelope and outflow cavity, suggesting dynamical or chemical interaction between outflow and 
envel ope. The effects o f outflows on the morphology of N2II"'" emission are also seen in other YSOs 



(e.g., IChen et al.ll2008l ). In addition, the emission extends to the very west of the observational 
field of view and becomes blended with noise. In this study we do not consider the western clumps 
(see Figure [D^b)); instead, we focus on the main part of the flattened envelope. 
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4. Analysis and Discussion 

4.1. The Column Density of the N2H+ Feature 

The N2H"'" molecule (diazenylium) has 15 J=l-0 rotational transitions; spectroscopically, seven 
sets of hyperfine com ponents are observe d because the lower state energy levels J=0 are degenerate 
within 10~^ H z (e.g., Daniel et al.ll2006l ). The relative rest frequencies of the components were de- 



termined bv lCaselli et al.l (|l995l ). The so-called isolated comp onent, the J.Fi ^ =101-012 component 
together with JFiF=101-010 and 101-011 at 93.176258 GHz (jhee et allboOll l. is not blended with 
the other sets of hyperfine components and is ideal for studying dynamics. We perform spectrum 
fitting for pixels with more than three S/N detection of N2H^ on the integrated intensity map, 
except the very west region close to the edge of the field of view, where the spectrum fitting fails to 
be reliable. For simplicity, we assume a uniform excitation temperature T^x = 10 K with the local 
thermodynamic equilibrium (LTE) approximation for the level populations. Also, Gaussian line s 
with the same line width are used (cf. GILDAS CLASS HFS procedure R and Ishirlev et"^l2005l). 



The isolated component accounts for a ratio of 3/27 of the total N2H"'" emission (jPaniel et al 



20061 ) , with which we convert our observed data to total emission accordingly. Then we calculate d 
the best fit values of opacity (r) (average r among seven hyperfine sets, see IShirley et al.ll2005l ). 
FWHM of velocity dispersion Aw, and fLSR foi' each pixel using the MATLAB function nlinfit%. A 
summary of the spectrum fitting as well as the uncertainty of fitting is in Table [2l In t he optically 



thin l imit and with the LTE ass umption, the column density can be estimated with (iMiao et al, 
199,4 iGoldsmith fc Lang^ll99d ^ 
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where 9a and 6b are observational beam size in arcsecond, v is the frequency in GHz, n is the dipole 
moment in debye, S is the line strength, is the upper state energy level, By(T) is the Planck 
function at temperature T, ly is the specific intensity in Jy beam~^, v is velocity in km s~^, Cr is 
the opacity correction factor 
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(jGoldsmith &: Langeiill999l ). Ej = J{J + l)hB is the energy level for the rotational transition, and 
B is the rotational constant. For linear molecule N2H"'", ^ is 3.40 debye and B is 46586.867 MHz 
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( JPL catalog: iPickett et alj 119981 : 1 Green et alj Il974l ) . Since the emission is more extended than 
the beam size, the beam fihing factor is assumed to be 1 across the map. Also, the background 
temperature Ti,g is assumed to be constant 2.73K. 

Altogether, the estimated N2H+ column density is shown by the dark contours in FigureEl with 
the velocity map shown in colormap in Figure El^b). The average column density is ~1.0xl0^^cm~^ 
with a peak of 
has a mean value of 



"^13qv-v-» 2 

■^3x1012 



3.8x10^"^ cm close to the protostar. The uncertainty varies across the map and 



cm 



The uncertainties of the column density estimation are dominated by the assumptions of spec- 
trum fitting. First, the largest error source for the derived column density is the uncertainty of the 
excitation temperature T^x- Typically, the excitation tem perature in dense cores and protostars 



Kirk et al. 


2007: 


Chen et al. 


2007) 



Benson et al.lll998l : IShirlev et al.ll200,4 



perature of L1157 was determined to be 8.9K bv lEmprechtinger et a.. 



(|2009| ). in general agreement 



with the rotational temperature of 13K obtained in lBachiller et al.l (|l993l ) using NH3 emission. In 
this study we assume a constant excitation temperature of lOK, and impose an uncertainty of 3K 
that is propagated with the other parameters. We used additional CARMA data that contain all 
hyperfine components in a lower spectral resolution for spectrum fitting with the CLASS package, 
and confirmed that the excitation temperature is around lOK near the protostar. 

Second, we assume that the N2H~'' lines are optically thin. In fact, the average opacity from 
our spectrum fitting is 0.31, consistent with the assumption of optically thin. Nonetheless, the 
correction factor helps mitigate the errors propagated to column density. In addition, any self- 
absorption effect is neglected given the optically thin assumption. However, the optically thin 
approximation becomes worse toward the center, where r=0.63ib0.03. 

LTE is assumed for the level populations of the rotational transitions. A constant proportion- 
ality between the hyperfine components is also adopted. The var iation of the r elativ e strengths 
between the hyperfine components as discussed in the appendix of IShirlev et al.l (|2005l ) is ignored 
in this study. 

For simplicity, we assume only one velocity component for the spectrum fitting, although this 
is not true in the inner envelope. We use a single Gaussian to fit the spectra over the map to 
estimate the N2H+ column density. The simplification is insufficient for the more complex velocity 
structures near the center, but the resulting difference is smaller than the uncertainty caused by 
the excitation temperature assumption. For the kinematics of the envelope, see later sections for 
further discussion. 

Unlike most N2H"'" studies, we fit the N2H^ spectrum using only the isolated component instead 
of all seven sets of transitions. The 2MIIz bandwidth in our observations cannot cover all seven sets 
simultaneously, so the excitation temperature distribution cannot be derived based on the relative 
strength between hyperfine components; therefore we assume a uniform excitation temperature. 
However, the difference should not be significant; it has been suggested that the two approaches 
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give results consistent within frequency resolution (e.g.. lMardones et al.lll997l : lEmpreclitinger et al 



20091 ). 



Indeed, our derived N2H"'" column density of L1157 is c omparable to other studies. For exa mple 



the column dens ity was reported to be 1.7x 10 cm in iBachiller &: Perez Gutierrea (|l997l ) and 



1.1x10 cm in lEmprechtinger et al.l (120091 ) . both using the IRAM 30m telescope with half-power 



beamwidth (HPBW) of 27" with different assumptions of excitation temperature. Our average 
column density (l.Ox lO^^cm^^) is slightly lower, as expected due to the interferometer resolving 
out large-scale emission. 



4.2. Correlation between N2H+ and 8 jj^m Absorption 



m 



The extended feat ure of N2H"'" emission coincides with the flattened absorption structure found 



Loonev et al.l (|2007l ). Figure [3] shows the absorption feature at 8fim overlaid with the contours 



of N2H"'" column density. At 8 fim, the polycyclic aromatic hydrocarbon (PAH) emission provides 
a bright background, against which the dust extinction is detected. Note that the spatial resolution 
of Spitzer IRAC (the diffraction limit is 1.8") is better than our N2H"'" observations (beam size is 
7.1"x6.3"). In general, the morphology of the absorbing dust is consistent with the N2H"'" gas, 
except the central region where the bipolar outflow feature dominates. From the 8 /im map, the 
ob served opacity is obt ained given the assumption of a constant background estimated off-source. 
In iLoonev et al.l (j2007l ). a simple model of the flattened envelope was constructed; however, the 
power-law index of the density profile could not be constrained. 



_The observed opacities were converted into total column density and mass in iLooney et al. 

= 5.912 cm^g-^ 



(120071) with the assumpt ion of a mass absorption coefficient of dust plus gas n%.o^r, 



fromlLi &: Draind ( 200ll ). The total absorbing mass was estimated to be 0.19 Mq in lLoonev et al 
(120071 ). For a comparison, the mass of the same extended region is '^0.76 Mq based on the derived 



N2II"'" column density and a hypothetical constant N2H+ abundance n(N2H+)/n(Il2) of 3.0x10^^*^. 
Note that the central part of the flattened envelope is excluded here. Although most concentrated, 
the central region is do minated by th e outfl ow activities at 8 /xm; therefore the method of optical 
depth determination in iLooney et al.l ( 2007 ) is not applicable for regions within ~8.4" of the pro- 
tostar. Also, the chemistry can become more complicated near the center so that the real N2II"'" 
distribution is not trivial. Furthermore, this comparison should not be overemphasized because of 
the following complexities: (1) an uncertain assumption of the N2II"'" abundance, (2) interferometric 
filtering, (3) the uncertainty of KH.Ojim,, and (4) ignorance of foreground emission when estimating 



mass using 8 fim absorption (jTobin et al.ll2009l ). 



Neverth eless, the tota l colu mn density, including gas and dust estimated from the 8.0/im 



extinction in ILooney et al.l (120071 ). can be compared with the N2II"'" column density derived from 
our spectrum fitting in this study. To facilitate the comparison for the flattened structure, we plot 
the profiles along the major axis of the flattened structure in Figure 4. Figure Sl^a) shows the 
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N2H"'" and estimated total column density, derived respectively from CARMA and Spitzer IRAC 
observations, for a cut of position angle 75° from n orth to east. The cir cles with error bars are the 



total column density from the absorption feature in lLoonev et alj (|2007l ). and the thick curves with 
shades are the N2H"'' column density from this study. 

As mentioned in the last section, the dominating uncertainty of the N2H"'" column density 
comes from the assumption of the excitation temperature Tf>x- In our spectrum fitting, a constant 
Tex of lOK is assumed. We estimate the uncertainty by imposing an uncertainty of excitation 
temperature ATea;=3K and the propagated error is shown as a color shaded region in Figure [H 
Because the column density is monotonically increasing with Tg^ (from Equations (1) and (3), 
A/tot oc Q{Tex)e''Tex cx Texe'^^ex )j the upper and lower bounds of the shaded region correspond to 
Te2;=13K and 7K, respectively. 



4.3. N2H+ Abundance 



We estimate the N2H"' " abundance by taki ng the ratio of [N2H"'"], from our CARMA observa- 



tions, and [H2], derived in iLoonev et al.l (j2007l ). along a cut of the extended envelope. In Figure 
mb), the ratio of [N2H+] to [H2] as a function of offset is plotted. Only the region away from the 
central outfiow-dominating zone and with more than three S/N of N2H"'" detection is considered 
here. While the radial profile of the 8 /xm absorption shows a good symmetry between east and 
west wings, the profile of N2H"'" does not. The N2H"'" column density drops in the west wing faster 
than the east. The average N2H"'" abundance is 3.0x10^^'' in the east wing, lower than the average 
abundanc e 5.7xl0~^'^ in the west wing . Our estimate of N2H"'" abundance is smaller than the value 
derived in lBachiller &: Perez Gutierrea ( 19971 . 3.8x10"^), while the difference can come from a dif- 



ferent estimate of H2 column density and interferometric filtering. Nonetheless, the overall N2H"'" 



prestellar cores and Clas s YSOs varies from 10 to 10 (e.g. 



Aikawa et al. 


2005: 


Maret et al. 


(e.g., 


Berein & Laneer 


1997). 



Enhancement of N2H^ abundance close to the center is seen in the east wing (Figure [D[|b)). 
This trend is also expected by chemical models. While the preced ing objects, starless cores, have 
been shown to have constant N2H"'" abundance (jTafalla et al.ll2002l ). the abundanc e profile evolves 
as a function of radius in the collapsing cores (e.g., the numerical chemical models: iLee et al.ll2004l : 



Aikawa et al.ll2005l : iTsamis et al.ll2008l ). In particular, the formation and de struction of NoH"*" is 



closely related to oth er species by the chemical reactions with N2 and CO (jWomack et al 



1992 



Jorgensen et al.l 120041 ). The dominant route to form N2H"'" is through 

H+ + N2^N2H++H2, 
while it is destroyed by dissociative electron recombination at high temperature 

N2H+ + e- ^ N2 + H, 



(4) 



(5) 
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and CO destroys N2H"'" in the gas phase by the reaction 



N2H+ + CO ^ N2 + HCO^ 



(6) 



Therefore, when C O starts to deplete onto dust grains at densities around 2-6 xlO^ cm ^ (e.g. 



Tafalla et al.ll2002l ) in the inner envelope, the abundance of N2H"'", as well as other nitrogen-bearing 
molecules, increases. 



On the other hand, depletion of N2H 
tion density is around 1 0^ cm~^ to 3x10^ cm 



(jBergin &: Langer 



1997 



is also expected in the inner envelope. The deple- 
^, above which the gaseous N2H"'" start to deplete 
Aikawa et al.l 120031 ) . However, whether depletion is seen in our observa- 



tions is not clear due to the fact that the region with density higher than the depletion density 
is not well resolved with the observational beam size. Also, information of extinction and total 
column density is missed owing to the bright outflow activity at 8 fim, making us unable to obtain 
the N2H"'" abundance near the center. Furthermore, the CO molecules, existing in the outflows or 
evaporated from the dust grains in the ver y inner envelope due to the heating from the central 
protostar, can destroy the N2H+ molecules (iLee et al.ll2004l ): the CO outflow and evaporation ef- 
fects can look similar to the depletio n effect froni the o bservational point of view. Eventually, a 
more careful study for L1157, such as lEvans et al.l (120051 ) for B335, a similar protostar at the same 
evolutionary stage as L1157, will be needed to consider multiple species and understand the system 
more thoroughly. 



4.4. Simple Dust Modeling 

While a large scale flattened envelope is detected by N2H+ emission, the dust emission at 3 
mm is compact and round; the goal in this section is to construct a model that can interpret both. 
To do so, we begin with fitting the dust continuum. We constructed a model that has a flattened 
geometry similar to the N2H"'" feature and predicts an observed spherical dust continuum. 

We have developed a radiative transfer code that considers density and temperature structures 
in three dimensions to do the dust continuum modeling, and compare with visibilities from inter- 
ferometric observations. First, a map of flux density is obtained with high numerical resolution 
given a model envelope. For each pixel on the pla.ne of sky, the flux is calculated by integrating the 



dust emission along the line of sight (e.g.. lAdamslll99ll ). The specific intensity can be expressed as 



f B^iT) e-^^ dr^ = f B^{T{r)) e'^^^"^ p{r) n^dr, (7) 

J los J los 

where B^{T) is the Planck function at dust temperature T, is the optical depth from the position 
r along the line of sight to the observer, p is the envelope density, and r denotes the position. T, 
p, and Tj, are all functions of r. In this study, the opacity of dust grains is assumed to be 
uniform throughout the envelope. We simulated all observational effects from the interferometric 
observations. The image is corrected by the primary beam effect according to the antenna size, and 
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Fourier transformed into visibilities with the observational u-v sampling. In the case of CARMA, 
the 6m BIMA dishes and 10m OVRO dishes give three different primary beams. Therefore we 
constructed separate primary-beam-corrected images for each kind of baseline, and sampled them 
with corresponding u-v spacings from real observations. However, the mosaic pattern was not 
considered in the dust modeling, although a five-pointing mosaic was used in the observations. The 
dust continuum is very compact with size much smaller than the FWHM of the primary beam, 
hence the effect is insignificant. Finally, the modeled visibilities are binned into u-v annuli around 
the source center and averaged vectorially. The flux density for each bin is compared with the 
observations, which are binned in the same way using the MIRIAD task UVAMP, as shown in 
Figure [5] by crosses with error bars. A Xr value is calculated to quantify the goodness of a fit. The 
comparison is done in visibility space so uncertainty from imaging processes such as the CLEAN 
algorithm is avoided. 

The code has been tested for accuracy and convergency. The results of spherica hy symmetric 



cases were checked by compar ing with those obtained by previous codes used in iLooney et al 



([2003) and IChiang et al.l (120081 ) . Because of the nature of higher density and higher derivatives in 
the inner envelope, smaller step sizes are required for the inner region. On the other hand, step 
size is not as sensitive for the outer region, so computing time is saved by using coarser steps. In 
the current code, three refinements of step sizes are used with extra fine steps taken for the very 
central region to achieve accuracy better than observational uncertainties. Also, cases of different 
numerical resolutions have been run to check convergency. In this study, the numerical resolution 
is 37.5AU, corresponding to 0.15" for L1157. Higher numerical resolution gives consistent fiux 
density within a difference smaller than 0.05 mJy. In addition, the code itself does not assume 
any symmetry. In other words, the model of envelope density and temperature can be arbitrary in 
three dimensions. 

With the observed fiattened structure in mind, we construct a model envelope of similar 
morphology. For simplicity, we adopted a power-law density and temperature profile. To impose 
the fiattened geometry, we use a spherically symmetric distribution and taper the density profile 
with a Gaussian mask diminishing along the minor axis of the disk orientation. In our model, 
the vertical scale height of the Gaussian is fixed to be 2000 AU for an ellipticity close to the 
observed fiattened feature. We orientate the fiattened envelope as the observed N2H"'" feature, 
that is, with a position an gle of 75°. We als o adopt an inclination angle of 80°, as determined 



from outfiow observations (jGueth et al.l Il996l ). However, the temperature distribution remains 
spherically symmetric. 

Note that this model is not physically motivated, but it interprets the observations with statis- 
tical significance. We adopt a power-law index of -2 for the dens ity structu re, similar to a singular 



isothermal sphere in the extreme case of a Bonnor-Ebert sphere (IShulll977l ). The dust temperature 



profile is a power law of index -2/(4-|-/3), where /? is the spectral index of opacity, from th e spher 



ically symmetric approximation assuming the heating is dominated by the central source (lAdams 



I991I ). We set the extrapolated temperature at lAU to be 300K and a lower limit of 5K representing 
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the external heating from cosmic ray or interstellar radiation. 



For the dust grain properties, a fiducial value o f the mass opacity coefficient n ^mm =0.0056 

0.01 cm2g-i fe.g.. lNatta et al.ll2004l ^ and (5 = 0.5 dKwon et al 



cm^g ^ is adopted, based on Kimm 



20091 ). The exact value of k is uncertain and can var y by an order of ma gnitude dependent on the 



gr ain properties (e.g., the inte rstellar MRN grains in lDraine &: Led 1 1984 and the coagulated grains 



m 



Ossenkopf &: Henninall994l ). Nevertheless, is expected to follow a power law with frequency v 
at long wavelengths, that is, k oc u^. The extrapolation at millimeter wavelengths is determined by 
the spectral i ndex /3, which is dependent on grain size, grain geometry, chemical composition, etc 
(jDraindl2006l ). For a young object like L1157, grain growth may have already started in the inner 
region. For simplicity, the dust grain properties is assumed to be uniform across the envelope and 
not depend ent on the radius. The possible radial gradient of grain properties, which was recently 
reported by lKwon et al.l (120091 ). is neglected in this study. 



We have explored two parameters: point-source flux (Fp), representing the flux contribution 
from the embedded circumstellar disk system inside the inner radius of the model envelope, and 
the extrapolated density at lAU (di), scaling with the total mass. A grid of model parameters 
is run with Fp from to 30 mJy in steps of 1 mJy, and di from 4x10^^^ to 20xl0~^^g cm~^ in 
steps of 0.5xl0~^^g cm^^. The inner and outer radius of the envelope are fixed to be 5AU and 
15,000AU, respectively. Varying the outer radius of the envelope does not make a big difference 
because little fiux is emitted from the low-density region compared to the denser region. Also, the 
interferometer resolves out some of the large-scale emission. The inner radius is correlated to the 
point-source flux and can be chosen as long as it is much smaller than the observational spatial 
resolution, while physically the dust destruction radius is of the order of O.IAU, inside which the 
temperature is too high for dust to exist. A summary of the Xr model fitting is shown in Figure 
[6l The filled circle marks the best fit and the contours show the parameter ranges with different 
confidence levels. The two parameters are correlated and cannot be clearly distinguished, because 
either a big density scaling factor di or a high point-source fiux Fp can result in a high peak of 
fiux density. The best-fit model has a point-source fiux Fp = 1 mJy and extrapolated density at 
lAU di = 1.35xl0^^^g cm^'^. A total envelope mass of 1.5 Mm is then implied in our model; this 
mass is general ly consistent with other mass estimates (e.g.. IShirley et al.ll2000l : iGueth et al.ll2003l : 



Froebrichll2005l ). The modeled flux density of the best fits as a function of uv distance is shown by 



the curve in Figure [5l 

We estimate the corresponding N2H"'" column density of the best-fit model obtained from the 
dust continuum fitting. To do the conversion, the dust-to-gas ratio is assumed to be uniformly 
1/100, as the typical value in the interstellar medium, and the N2H"'" abundance is assumed to be 
uniformly 3.0x 10~^*^. Further, we simulate the depletion effect by introducing a threshold depletion 
density, above which the N2H"'" density does not increase with H2 density. We adopt a depletion 
density of 1.5xl0^cm~^. To compare with observations, N2H^ volume density is integrated along 
the line of sight to calculate the column density across the map. The model N2H"'" column density 
and the detected N2H"'" emission are shown in Figure [71 
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As seen in Figure [71 our best-fit dust model gives a consistent map of N2H"'" column density. 
This simple model illustrates how we do not observe extended dust continuum while the N2H"'" 
emission is detected to be extended more than 80" across. The dust emission, which depends on 
both density and temperature, in the outer region of the flattened envelope is too weak to be 
detectable by these observations. 



4.5. Gas Kinematics 

The complex N2H^" spectral information reveals a composite system of multiple dynamic com- 
ponents (Figure [3^b)). To understand the detailed kinematics, we perform a simple analysis using 
position-velocity (PV) diagrams. Figure [8|^b) shows the PV diagram for the N2H"'" JFiF=Wl-012 
component along the major axis of the flattened envelope at position angle 75° from north to east. 
This transition is isolated so no confusion is caused due to blending with other hyperfine lines. 

We construct a simple model to simulate the PV diagram assuming various velocity structures 
of the envelope. Note that the kinematics modeling is independent of the dust modeling, except 
that the best-fi t density structure fr om the previous section is used. We adopt the procedure 



as discussed in lOhashi et al.l (|l997al ) and consider a spatially thin cut along the major axis of 
the envelope. Both the inclination and opacity effects are ignored. The column density in each 
velocity channel as a function of offset is estimated by integrating density along the line of sight and 
sorting into observational spectral bins by the projected velocity. The depletion effect is simulated 
by assuming constant N2H"'" density for the central region where H2 gas density is higher than 
1.5xl0^cm~^, the same depletion density as adopted in the previous section. With the depletion 
effect, the actual density structure of the inner envelope does not play a large role. The spatial 
distribution is then convolved with the observational beam. However, unlike the dust modeling, 
the interferometric filtering effects are not taken into account for the kinematics study. Nonetheless 
the large-scale emission from background clouds is not included; this is similar to interferometry 
resolving out large scale emission. Finally, the modeled results are shown by contours in PV 
diagrams and compared with observations. 



.5.1. Spectroscopic Signatures of Infall Motion 



The double-peaked feature seen in the central region indicates multiple velocity components 
in the protostellar envelope (Figure Ell^a)). The redshifted and blueshifted peaks have a velocity 
difference of ~0.4 km s^^, inferring a relative motion between two parts of the envelope. One 
interpretation for this is the infall motion of the inner envelope, which can also be the cause of 
the high- velocity wings. Gr avitationa l collapse takes place in the early stage o f star formation 



(e.g., theoretical studies in IShu 



1977; 



Hunteii I1977I : iTassis &: Mouschoviad l2005l ) . However, the 



predicted infall velocity structures are very different from model to model, especially in the inner 
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envelope. Detailed comparison of various theoretical models with observations is beyond the scope 
of this paper and may be studied with observations of higher spatial and spectral resolutions in 
the future. Nevertheless, infall is expected in the in ner envelope of L1157. Previously, L1157 has 
been identified as an infall candi date by single dish (IGregersen et al.l 119971 : iMardones et al.l 119971 ) 
and interferometric observations (jGueth et al.lll997l : iBeltran et al.ll2004l ). In addition, detection of 
methanol from the accretion shocks on the sr nall scale also support s the picture of envelope material 
infalling onto the forming protostellar disk (jVelusamy et al.l 12002 ). 



Figure [9)[a) shows the modeled PV contours for a pure infall motion. A simple velocity profile 
V ~ r~^/^, representing ideal free-fall, is used. The simpl e infall model shows a double-peaked 
feature toward the center of the envelope, as obtained by lOhashi et al.l (|l997al ). However, the 
model cannot explain the asymmetry of the peaks. A more sophistic ated model tha t cons iders full 
geometry, radiative transfer, and different theoretical profiles such as Momose et al. (jl998l ) may be 
necessary. 

We exclude this double-peaked feature as arising from self-ab sorption, commonly used as an 
indicator of infall motion in the prestellar cores (e.g., lEvanslll999l ). In such an approach of iden- 
tifying infall candidates, a stronger blueshifted peak and a weaker redshifted peak are expected 
for infalling sources observ ed by an opaque l i ne because the front half of the cloud causes a red- 



shifted absorption dip (e.g.. lMyers et al.lll996l : iMasunaga &: Inutsukall2000l ). while an optically thin 
line should be observed within two velocity peaks. In this specific case, since the N2H"*" lines are 
moderately optically thin, this scenario of self-absorption is ruled out, although a similar profile is 
observed. However, the spectrum is sensitive to the detailed envelope model such as optical depth, 
turb ulent yelocity dispersion, spatial structures of ve locity, density, and excitation temperature, etc 



(e.g., lzhoulll992l : IWard-ThomDson fc BucklevlboOlh . The complexity makes it difficult to justify 



the real physical properties causing the observed spectra. In addition, the interferometry selectively 
observes the small-scale emission, leaving out the larg e-scale stati c cloud in the outer envelope and 
resulting in spectrum with a deeper zero- velocity dip (|Choill2002l ). 



While a double peaked profile of N2H~'" is seen in our interferometric observations. IMardones et al 



(jl997l ) observed the same line using IRAM 30m telescope, and found it to be single-peaked and 
used it as the optically thin reference to compare with other optically t hick lines for studying infall . 
The same source was observed again by the IRAM 30m telescope in lEmprechtinger et al.l (|2009l ) 
and similar results were obtained. [§ These support our spectrum fitting results that the N2H"'" 
J=l-0 lines are optically thin. The major differences of our observation are a smaller beam size and 
the aperture synthesis. Therefore the cause of discrepancy is deduced to be (1) a beam-smearing 
effect and/or (2) interferometric filtering. While the overall spectrum is dominated by the inner en- 
velope, the small-scale structures are not distinguishable by the single-dish observations. However, 
observations of higher resolutions with interferometers reveal other complexities. For example, the 



^However, our recent spectral data taken by IRAM 30m seems to show a double-peaked feature in the strongest 
hyperfine component (J.J. Tobin et al. in preparation). More careful examination is in process. 
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degre e of profile asymmetry can be enhanced if observed with a higher angular resolution (jChoi 
2002 ) . Also, we test our modeling routine by convolving the spectrum with different beam sizes, 
and confirm the beam-smearing effect. A double-peaked feature can look like a single Gaussian if 
the observational beam size is not able to resolve the infall radius. 

Second and more importantly, interferometric observations reveal only the targeted protostel- 
lar envelope at the expense of missing flux f rom large-sc ale structures, arguably dominated by 
the static foreground and background clouds. IChoil (j2002l ) has shown that missing short spacing 
flux can affect the self-absorption dip for optically thick lines. It can cause similar effects for 
optically thin lines too. In particular, large scale background material is mostly static compared 
to the infalling inner envelope and contributes more flux at the LSR velocity. Nonetheless, the 
material in the large scale cloud does not participate in star formation as actively as the inner en- 



velope. For comparison with the spectrum of IRAM 30m single-dis h observations (iMardones et al 



Bachiller &: Perez Gutierrez 



1997 



19971 

27" Gaussian beam (the beam size of IRAM 30m observations) and found that 



Emprechtinger et al.ll2009l ). we smooth our images with a 

20% of the total 



flux is filtered out by CARMA, also causing a dip around the LSR velocity. 

A large line width of N2H+ J=l-0 is found by the single dish observations. The fitted line 
widths of the single-p eaked spectra are 0.71 and 0.65 km s"^ for the IRAM 30m observations in 
Mardones et al.l (jl997l ) and lEmprechtinger et al.l ([20091) , respectively. We compare the single-dish 
observations with our smoothed spectrum. If we fit our double-peaked spectrum with a single 
Gaussian, the FWHM is 0.88±0.05 km s'^ But if we fit the spectrum with two Gaussians of the 
same width, the FWHM is 0.44 ±0.02 km s~^, comparable to the sonic line width ~0.45 km s~^ 
and larger than the thermal line width ~0.13 km s^^ with the assumed temperature of lOK. The 
dynamics at small scales can be obscured in single-dish observations but can be revealed by high- 
resolution interferometric observations. If this phenomenon is common to other sources, it may be 



why the average line width of Class YSOs is la rge (~0.61 km s ^, Emprechtinger et al.ll2009l ) 



compared to that of starless cores (0.2-0.4 km s~^ lLee et allboOlh . 



We cannot exclude the possibility that the double-peaked feature is caused by factors other 
than infall motion in the inner envelope. For example, outflows contaminate the mol ecular tracers 



of env elope material. A similar case of the Class protostar B335 was studied by IWilner et al 



(120001 ) . While the single-dish observation matched well with the inside-out infall model (jZhou et al 



19931 ). the interferometric observation brought up a more realistic scenario. The optically thin CS 
J =5-4 line is shown to be dominated by small-scale outflow clumps in high-resolution observations 
(jWilner et al.l I2OOOI ) . A single-peaked spectrum seen in low-resolution observations can contain 
two velocity components and hence a double-peaked feature is shown with interferometric filtering. 
However, some molecular lines do not trace outflows as closely as other species. It has been 
sugge sted that N2H"'" lines trace the quiescent cores but not the shocked outflow gas (IBachiller 
19961 ). The CO molecules in the outflows can destroy the N2H"'" molecules (Equation (6)). Indeed, 
while the abundan ce of some species is enhanced by the outflow shocks, N2H"'" is not detected at 
the shock regions (jBachiller &: Perez Gutierrez! 119971 ). In other words, the outflow contamination 
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is minimized for this specie, although some outflow effects are unavoidable (for example, see the 
outflow-envelope interface at the southeast extension in Figure [T]^b)). 

The double-peaked spectrum can also be caused by unrelated dense clumps that happen to be 
in the same line of sight. This is less likely because channel maps do not show traits of unrelated 
components. If clumps exist in the foreground or background, they are more likely to show up 
individually at a peculiar velocity. In contrast, the morphology of the envelope at different velocities 
are systematically consistent in our observations. Also, N2H"'' does not pick up other velocity 
components as easily as CO because it requires high densities to form. Infall motion, not necessarily 
gravitational infall, is suggested. 



4-5.2. Rotation 



While axisymmetric infall cannot explain the velocity difference between the east and west 
extensions of the envelope, the differential velocity may come from rotation of the envelope. In this 
section, we test some simple rotation curves. 

First, we test Keplerian rotation. If the dynamics of the envelope is dominated by Keplerian 
rotation around a central point mass much larger than the cumulative envelope mass, the central 
mass has to be smaller than 0.1 Mq to explain the observed velocity differential of the envelope. 
In Figure E^b), the pure Keplerian curves for a point mass of 0.01, 0.02, and 0.04 are shown. 
Apparently, the deduced central mass is around or even smaller than the mass of the envelope, 
which means that the cumulative envelope mass cannot be neglected. This is expected because for 
Class YSOs, ma ss is mostly distributed in the envelope rather than the central protostar (e.g., 
Andre et al.lll993l ). Next, we consider complete Keplerian rotation with cumulative envelope mass. 
The velocity is {GMr/R)^/'^, where Mr is the mass contained within the radius R. An envelope with 
p oc gives a flat rotation curve. The best-fit density obtained in the previous dust continuum 
fitting gives a roughly constant rotation velocity 0.5 km s^^, much larger than the observed value. 
These exercises show that Keplerian rotation is ruled out for the large-scale envelope. 

In fact, the system is too young to construct large-scale Keplerian rotation. For example, the 
period of rotation with velocity 0.2 km/s at 1000 AU is --^1.5x10^ yr, larger than the typical age 
of Class protostars. As seen in the PV diagram, the size scale of rotation is around 10,000AU 
and requires even longer dynamical timescale. Again, the dynamical time estimation suggests the 
unsuitability of pure Keplerian rotation for this system. 

Another possible scenario is that t here exists solid-body rotation in L1157. This has been seen 

The solid-body rotation is probably 



in other Class YSOs such as HH212 (Wiseman et al 



2001 



inherited from the initial conditions of the large-scal e clouds or filaments . Angular momentum plays 
an important role in the protostar evolution (e.g., iBodenheimeri Il995l ). In particular , the initial 



condition of rotation is closely related to the core morphology and fragmentation (e.g., ISaigo et al 



20081 ). The extended envelope of L1157 has a velocity gradient of around 1.5 km s ^ pc ^ assuming 
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solid-body rotation (shown by the dashed hne in Figure [Hl^b) and (c)). The velocity gradient is 



much smaller than what was found in the survey of Class protostars (~ 7 km s ^ pc ^. IChen et al 



2007); instead, it resernbles the typical ve locity gradient found in the dense clouds (1-2 km s ^ pc 



Goodman et al.lll993l : ICaselli et al.ll2002l ). Moreover, the rotating N2H"'" envelope of L1157 has a 
size scale much larger than the typical size scale of collapsing envelopes. Both the large size and 
the slow bulk rotation imply properties more similar to prestellar cores than collapsing envelopes. 
After the protostar has formed in the central densest region, the kinematics of the envelope can 
still be dominated by the remnant rotation of the parent dense cloud at large scales, while the infall 
motion takes over at small scales. The good alignment of the rotation with other features such as 
flattened geometry and outflow direction also suggest a consistent picture. 

Considering the envelope with radius of 10,000AU and the fitted solid-body rotation, the 
specific angular momentum is around 4xl0~^ km s~^ pc. We can locate it at a specific angu- 
lar momentum-rotation radius plot to compare with the other dense cores and protostars (e.g., 
Ohashi et al.l Il997bl : iBelloche et al.l I2OO2I : IChen et al.l 120071 ) . Although as an infalling envelope at 
small scales, the large-scale properties of the L1157 flattened envelope resemble those of prestellar 
cores more than those of more evolved protostars. Additionally, the fitted angular velocity implies 
that if the large-scale envelope collapses to form a circumstellar disk, the envelope material would 
fall onto the midplane within a centrifugal radius of Tc = R^Q? /GM ~ 500 AU. This radius is much 
larger than the observed T Tauri disks. All above implies that the large-scale flattened envelope is 
probably not involved in the dynamical infall activity. 

The interpretation of solid-body rotation is well suited from many theoretical points of v iew. 
For example, magnetic braking can induce solid-body rotation (IBasu fc Mouschoviaa 11994). A 
turbu lent core can also yield velocity gradient similar to uniform rotation (IBurkert &: Bodenheimer 
2000). Here we calculated the Prot parameter, defined as the ratio of rotational kinetic energy to 



gravitational energy (3, 



•rot 



7T-, 



qGM'^/R ~ 2gGM ' different radii of the flattened envelope. (Note 
that this Prot is different from the opacity spectral index mentioned in the previous sections.) For 
L1157, Prot is smaller than 2% throughout the envelope, suggesting that the flattened structure is 
not supported by rotation. 

The combined best-fit model of simple infall and solid-body rotation is shown by the black 
contours in Figure [9l^c) . Both the double-peaked feature and the large-scale velocity gradient can 
be explained by this model. However, the asymmetric features are impossible to model with an 
axisymmetric model in the optically thin case. Local dumpiness can be the cause of asymmetry. 
Without further observations with higher spatial and spectral resolutions it is hard to justify the 
detailed properties. 
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4.6. Overall Gas-Dust Comparison and Global Picture of L1157 



By considering the dust and gas information of Class YSO LI 157 altogether, we are able 
to construct an overall physical picture of the system. The dust absorption at 8 /im. from Spitzer 
shows a large-scale extended dense cloud that is flattened perpendicular to the outflow direction, 
while the 3 mm dust emission from CARMA D- and E-array shows a compact spherical structure. 
These two observations of dust detect different components of the envelope. The extended envelope 
detected by 8 /xm absorption is likely part of the parent cloud, while the very inner region seen by 
our dust continuum is the collapsing envelope. The flattened geometry may result from physical 
processes of core formation. On the other hand, the N2H"'" emission provides a consistent view both 
morphologically and kinematically. The N2H"'" feature coincides with the extended dust absorption 
seen at 8 fim, meaning that the same cloud component is observed. A slow solid-body rotation at 
large scales is seen along the major axis of elongation, but the flattened structure is not supported 
by rotation. This parent cloud resembles the physical properties of a prestellar core, while the 
innermost region is decoupled and undergoes (gravitational) collapse. 

A simila r dynamical scenari o has also been suggested for another Class object, IRAM 
04191+1522 (IBelloche et ahlbood ). They have obtained a decoupling radius of ~3500AU, which 
divides the envelope into a rapidly rotating inner region with free-fall motion and a transition region 
connected to the ambient slowly rotating cloud. Their estimated angular velocity at large scales 
(~1.9 km s~^ pc~^ at 7000 AU and <0.5-l km s~^ pc^^ at 11,000 AU) are comparable to what we 
find for L1157. More recent interferometric observations revealed faster rotation at smaller scales, 



but resolved out the large-scale structures (jBelloche &: Andrei 12004 : iLee et al.ll2005l ). 



We construct an edge-on disk-like envelope model that fits our L1157 data. Given the estimated 
column density and assumptions of envelope geometry, the volume density averaged along the line 
of sight can be estimated. With the assumption of an outer radius of 15,000AU and a constant 
N2H+ abundance of S.Ox lO"^'', as used in our dust modeling, the average gas volume density ranges 
from ~2xl0^ to '^3xl0^cm~^. This volume density is consistent with the detectable density of 
N2H"'", that is, it is close to the critical density lO^cm"'^. On the other hand, the average density 
should not be higher than the depletion density (~10®-10^ cm~^); this can give a constraint on the 
thickness of the cloud. 

The geometric structures of L1157 at different size scales are shown at various observations 
(Figure [T0|) . At s mall scales, L1157 shows a nearly sphe rical morphology, as commonly seen in 
other YSOs (e.g., iLoonev et al.l I2OOOI : IShirlev et al.ll2000l ): detailed comparisons between obser- 
vations and theorie s of collapsing envelopes are usually done assuming spherical symmetry (e.g., 
Looney et al.ll2003l ). In particular, high-resolution observations done by interferometry only show 
the contributions from the inner envelope while the large-scale structures are mostly resolved out. 
Deviation from spherical symmetry can be significant, especially at large scales. At a younger 



stage of evolution, large-scale dense cores are commonly observed to be elongated (e.g.. lMyers et al 



1991 



). Although the underlying dynamical processes of the flattened cores are controversial (e.g. 
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Gammie et alj l2003l : iTassisI 120071 : lOffner &: Krumhola 12003 ) , n on-spherical structure in the initial 



cond ition can play an important role for their evolution (e.g., iGalli &: Shulll993l : iHartmann et al 
19961 ). For L1157, both the N2H"'" feature and 8 fim absorption reveal the flattened structure at 
the size scale of ~0.1 pc, comparable to a typical prestellar core. At a larger scale, the morphol- 
ogy becomes even more irregular. For example, the Digitized Sky Survey (DSS) |I| optical image 
of L1157 shows a large-scale structure of irregular shape. This is more like the initial condition 
but not the consequence of star formation. The cloud is filament-like and irregular at very large 
scale, while a nearly spherical inner envelope is embedded in the intermediate-size, flattened outer 
envelope extended perpendicular to the outflows. 

L1157 is an interesting and relatively simple case for early star formation. While many Class 
YSOs are associated with large-sc ale non-axisymme tric outer envelopes, the geometry of L1157 
is highly flattened and symmetric (iTobin et al.M2009l ). L1157 is a perfect site for observational 
experiments, presenting a typical Class YSO with less complexity. The properties of L1157 
may be generalized to other Class YSOs. For example, the Class collapsing envelope with its 
embedded protostar is surrounded by the outer envelope left by its parent core that may or may 
not gravitationally collapse at a later time. In the case of L1157, we are able to detect this outer 
envelope in a flattened structure at both 8 /um and N2H"'" emission, revealing a phase that shows 
both properties of the preceding prestellar core stage at large scale and properties of the current 
Class stage at small scale. In other cases, the outer envelope can become more complicated that 
a non-axisymmetric model will be needed to interpret the observations. 



5. Summary 

1. We observed the dust continuum and N2II"*" gas emission at 3 mm toward the Class YSO 
L1157 IRS with an angular resolution of ~7" using CARMA at D- and E-array. Spectra of 
the N2II"'" isolated component JFiF=101-012 were obtained with resolution of ~0.1 km s~^. 
While the 3 mm dust continuum detects the inner envelope, which is compact and nearly 
spherical, the N2H"'" emission shows a huge flattened structure with a lir iear size of ^2 ,000A U, 



coinciding with the disk-like feature found by the 8 fim absorption in lLooney et al.l ( 20071 ) 



2. By fitting the spectra, we estimated the gas column density across the flattened envelope 
and compared with the dust column density deduced from the 8 /im absorption feature. We 
derived the N2H"'" abundance and found results consistent with what is expected from chemical 
models. Further, we examined the radial profile of N2H"'" abundance along the major axis of 
the fiattened envelope, and showed the asymmetry between the east and west extension. 



^The Digitized Sky Survey was produced at the Space Telescope Science Institute under U.S. Government grant 
NAG W-2166. The images of these surveys are based on photographic data obtained using the Oschin Schmidt Tele- 
scope on Palomar Mountain and the UK Schmidt Telescope. The plates were processed into the present compressed 
digital form with the permission of these institutions. 
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3. We constructed a simple flattened envelope model that fits the compact dust continuum; fur- 
ther, the model-expected gas column density is consistent with the extended N2H"'" emission. 
It follows that the deviation from spherical symmetry can be important at large scales for 
protostellar envelopes. However, this model is not motivated by a theoretical model. 

4. We did PV contour modeling and studied the kinematics of the N2H+ feature. The spectrum 
of the central part of the system shows a double-peaked feature, implying infall. The large- 
scale component can be described by slow solid-body rotation comparable to the properties 
of a typical prestellar core. This large-scale filament may arguably result from the dynamical 
processes in the early core formation, while only the very inner part is actively involved in 
the protostar formation. 
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Fig. 1. — (a) 3 mm dust continuum of L1157 (contours and gray scale). The contour levels 
are [-2, 3, 6, 10, 14, 20, 26] xo", where a=l mJy beam~^ is the noise level, and the beam size is 
7.3" X 6.5" at a position angle of 86.1°. (b) The integrated intensity of N2H+ over the isolated 
hyperfine component JFiF=Wl-012 (2.21-3.58 km s~^, contours and gray scale). The star marks 
the position of the central protostar. The contour levels are [3, 6, 10, 15, 20, 30, 40, 60, 80]xc7, 
where (T=0.01 Jy beam~^ km s^^, and the beam size is 7.1" x6. 3" at a position angle of 86.6°. 
Although negative contours exist due to the filtered-out large-scale structures, they are not plotted 
for simplicity. 
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Fig. 2. — Channel maps of L1157 for the N2H^ Ji<'ii<'=101-012 transition (contours and gray 
scale). The numbers in the upper-right corner are the LSR velocity for each channel. The star 
marks the position of the central protostar. The contour levels are [2, 3, 4, 5, 6, 7, 9, 11] xcr, where 
a =0.08 Jy beam~^, and the beam size is 7.1"x6.3" at a position angle of 86.6°. Negative contours 
are not plotted for simplicity. 
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Fig. 3. — (a) IRAC 8 fim image of L1157 (iLooney et al.ll2007l . grayscale and gray contours) overlaid 
with N2H"'' column density {red contours) derived from the CARMA observations. The star marks 
the position of the protostar. The {red) contour levels for the N2H^ column density are [0.1, 0.5, 
1.0, 1.5, 2.0, 2.5, 3.0, 3.5] X lO^'^cm"'^ (see the text for discussion of uncertainty); and the {gray) 
contours showing the absorption features are [7, 8, 9, 10] x a, where a = 0.042 MJy~^sr~^ is the 
noise level. The circle at the lower-left corner shows the beam size of the CARMA observations, 
(b) Velocity map (vlsr) overlaid by the column density contours from (a), both derived from the 
CARMA observations of the N2H+ flattened envelope. 
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Fig. 4. — (a) Column density of N2H"'" (from our CARMA observation, curve with shading) and 
total gas plus dust (from Spitzer observation, circles with error bars) along the major axis of the 
extended envelope, (b) N2H+ abundance profile by taking the ratio of [N2H+] and [H2] from the 
CARMA observations and 8 /im absorption feature. The shaded region shows the uncertainty. 
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Fig. 5. — Annuli-averaged flux density for the 3mm dust continuum of L1157 {crosses with error 
bars). The curve show the best-fit model binned in the same fashion. The dust model has a power- 
law density profile with index = -2 tapered by a vertical Gaussian with a scale height of 2000AU. 
The extrapolated density at lAU is 1.35x10"^^ g cm~^ and an unresolved point source fiux of 
ImJy is added. The inner radius of the envelope is 5AU, and the outer radius is 15000 AU. The 
reduced Xr is 0.225 for this best-fit model. 
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Fig. 6. — Summary of the dust continuum fits. Circles show the grid of models with sizes and 
shades representing the goodness of fits. Contours are confidence levels of 99%, 95%, 90%, 80%, 
and 50%. The filled circle at 1.35xl0~^^g cm~^ and 1 mJy is the best fit. 
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Fig. 7. — N2H+ column density of the best-fit model by dust continuum modeling {black con- 
tours) and observations {grey contours). The contour levels (for both) are [0.1, 0.5, 1.0, 1.5, 2.5, 
3.5] X lO^^cm^^. The depletion effect is taken into account at the central region where the density 
is above 1.5xl0^cm~^. 
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Fig. 8. — (a) Spectrum of the N2H+ JFiF=101-012 line at the location of the protostar. (b) 
Position-velocity diagram of the N2H+ JFiF=101-012 line along the major axis of the flattened 
envelope feature, a slice with position angle of 75°. The contour levels are [2, 3, 4, 5, 6, 7, 9, 11] xcr, 
where a =0.08 Jy beam~^. The angular and spectral resolutions of observations are shown, and 
the pixel size of the gray scale is the interferometric imaging cell size. 
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Fig. 9. — Observed i^gray contours) and model {blue contours) PV diagrams for the flattened 
envelope of L1157: (a) pme infall motion, (b) solid curves show the point-mass Keplerian rotation 
for M^f=0.01, 0.02, and 0.04 Mq, dashed line shows solid-body rotation with angular velocity of 1.5 
km s~^ pc~^ , and (c) combined model of infall plus solid-body rotation. 
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Fig. 10. — L1157 at different scales seen by various observations: (a) DSS optical image, (b) Spitzer 
IRAC 8 image, (c) zoomed-in 8 fim image overlaid with the N2H+ column density map, (d) 
3 mm dust continuum observed by CARMA. Different geometric characteristics are shown. The 
cloud is irregular at scale of ~0.5 pc, flattened perpendicular to the outflows at scale of ~15,000 
AU, and nearly spherical at size scales smaller than ~5,000 AU. 
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Table 1. Summary of 3 mm Dust Continuum 



Parameter Value 

Peak flux density (mjy beam"^) 30.3 

Noise level <j (mJy beam~^) 1.0 

Total flux density (mJy) 48.3 ± 0.7 

Gaussian fit 

Major axis (arcscc): 9.4 ± 0.9 

Minor axis (arcscc): 8.1 it 0.4 

Position angle (degree) -68 it 18 

Deconvolved Size 

Major axis (arcsec) 6.1 it 2.3 

Minor axis (arcsec) 4.4 it 0.8 

Position angle (degree) -54 it 21 



Note. — The results of Gaussian fit with the 
dust continuum done by the MIRIAD task IMFIT. 



Table 2. Summary of N2H+ Spectrum Fitting 



Parameter Gentral Pixel Mean Value Mean a 

fLSR (km s-l) 2.85 2.82 0.02 

A?;(km s'^) 0.80 0.44 0.03 

< T > 0.63 0.31 0.03 

Derived values : 

Column density(10i2 cm-2) 38.3 9.7 3.2 



